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Abstract. The relationship between the properties of a dynamical system 
and the structure of its defining equations has long been studied in many 
contexts. Here we study this problem for the class of conjunctive (resp. dis- 
junctive) Boolean networks, that is, Boolean networks in which all Boolean 
functions are constructed with the AND (resp. OR) operator only. The main 
results of this paper describe network dynamics in terms of the structure of 
the network dependency graph (topology). For a given such network, all pos- 
sible limit cycle lengths are computed and lower and upper bounds for the 
number of cycles of each length are given. In particulax, the exact number 
of fixed points is obtained. The bounds are in terms of structural features of 
the dependency graph and its partially ordered set of strongly connected com- 
ponents. For networks with strongly connected dependency graph, the exact 
cycle structure is computed. 



1. Introduction 

The understanding of the relationship between structural features of dynamical 
systems and the resulting dynamics is an important problem that has been studied 
extensively in the dynamical systems literature. For example, work by Golubitsky 
and Stewart [l[ about coupled cell dynamical systems given by coupled systems of 
ODEs attempts to obtain information about system dynamics from the topology of 
the graph of couplings. Albert and Othemer [4] used a Boolean network model to 
show that the expression patterns of the segment polarity in Drosophila are deter- 
mined by the topology of its gene regulatory network. Thomas et al 0] conjectured 
that negative feedback loops are necessary for periodic dynamics whereas positive 
feedback loops are necessary for multistationarity. These conjectures have been the 
subject of many published articles [3;[l;[ll. In [3] we demonstrated that networks 
with a large number of independent negative feedback loops tend to have few limit 
cycles and these cycles are usually long. 

In this paper we study the effect of the network topology on the dynamics of a 
family of Boolean networks. Boolean networks in general, and cellular automata in 
particular, have long been used to model and simulate a wide range of phenomena, 
from logic circuits in engineering and gene regulatory networks in molecular bio logy 
0; H; (E [l^; [HI; 12 1 to population dynamics and the spread of epidemics fl3: 14|. 



Especially for large networks, e.g., many agent-based simulations, it becomes in- 
feasible to simulate the system extensively in order to obtain information about 
its dynamic properties, even if such simulation is possible. In such instances it 
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becomes important from a practical point of view to be able to derive information 
about network dynamics from structural information. But the problem is of interest 
in its own right, in particular in the more general context of time-discrete dynami- 
cal systems over general finite fields. These have been studied extensively, see, for 
examp^leJllJliJl^^ . They have a wide range of applications in engineering 



13; [iO; [iliniiliiniiniHand recently in computational biology [H; H 

Let F2 := {0, 1} be the Galois field with two elements. We view a Boolean 
network on n variables as a dynamical system 

Here, each of the coordinate functions fi : — > F2 is a Boolean function which 
can be written uniquely as a polynomial where the exponent of each variable in each 
term is at most one |28|. In particular AND{a, b) = ab and OR{a, b) = a + b + ab. 
We use polynomial forms of Boolean functions throughout this paper. 

Two directed graphs are usually assigned to each such system: The dependency 
graph which encodes the static relationships among the nodes of the network and 
the phase space which describe the dynamic behavior of the network. In this paper 
we focus on the question of deriving information about the phase space of a Boolean 
network from the structure of its dependency graph. Next we define these graphs. 

The dependency graph of / has n vertices corresponding to the Boolean 

variables of /. There is a directed edge i — > j if Xi appears in the 

function fj. That is. T){f) encodes the variable dependencies in /. It is similar to 
the coupling graph of Stewart and Golubitsky. 

Example 1.1. The dependency graph of / = {x2X3,xi,X2,X3X4,XiXg,X3X4X5) : 
¥2 — > F2 is the directed graph in Figure [TJ 




Figure 1 . The dependency graph of / from Example 11.11 

The dynamics of / is encoded by its phase space, denoted by S{f). It is the 
directed graph with vertex set Fj and a directed edge from u to v if /(u) = v. For 
each u G Fj , the sequence {u, /(u), ./^(u), . . . } is called the orbit of u. If u = ./*(u) 
and t is the smallest such number, the sequence {u, /(u), /^(u), . . . , /*^^(u)} is 
called a limit cycle of length t denoted by Ct, and u is called a periodic point of 
period t. The point u is called a fixed point if /(u) = u. If every limit cycle is of 
length 1, the system / is called a fixed-point system. Since F2 is finite, every orbit 
must include a limit cycle. We denote the cycle structure of / in the form of the 
generating function 

00 

(1.1) C(/) = ^C(/).C., 
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where C{f)i denotes the number of cycles Ci of length i in the phase space of /. 
Since the phase space of / is finite, C{f)i = for almost all i. 

The height of /, denoted by h(/), is the least positive integer s such that /''(u) 
is a periodic point for all u G . A component of the phase space S{f) consists of 
a limit cycle and all orbits of / that contain it. Hence, the phase space is a disjoint 
union of components. We define the period of /, denoted by p(/), to be the least 
common multiple of the lengths of all limit cycles in the phase space of /. 

Example 1.2. Let / : — > Fj be given by f{xi,X2,X3) = {x2X3,xi + x^, xiX2)- 
The phase space of / has two components, containing one fixed point and one limit 
cycle of length two, that is C{f) ~ Ci + C2, see Figure[2l It is clear from the phase 
space that h(/) = p(/) = 2. The phase space in Figure[2]was generated using DVD 




Without exhaustive iteration and just by analyzing what can we say about 
Namely, what is the period of /, the height of /, or the generating function 
C(/)? This question is NP-hard in general, so it is important to limit the class of 
Boolean networks considered. Next we list some of the main known results. 

• When all coordinate Boolean functions are the XOR function (that is, the 
functions arc linear polynomials), the above questions have been answered 
completely. In fact the questions have been answered for linear systems 
over any Galois field [13; 0; [3^. We have developed and implemented 
algorithms that answer the questions above, see 31|. 

• For Boolean networks where all coordinate functions are symmetric thresh- 
old functions, it has been shown that all cycles in the phase space are either 
fixed points or of length two js^l , 

• For Boolean cellular automata with the majority rule, the number of fixed 
points was determined in [3^, and 

• In [s^l, the authors studied AND-OR networks (Boolean network where 
each local function is cither an OR or AND function) with directed de- 
pendency graphs. Formulae for the maximum number of fixed points are 
obtained, and 

• The main result in [s^t is an upper bound for the number of fixed points 
in Boolean regulatory networks. 

In this paper we focus on the class of conjunctive and disjunctive Boolean net- 
works, that is, Boolean networks where all of their coordinate functions are either 
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the AND function or the OR function. The foUowing represent main previous 
attempts to mathematically analyze this class of Boolean networks. 

• A family of this class of networks have been analyzed in [s^. The au- 
thors studied conjunctive Boolean networks (which they called OR-nets) 
and disjunctive Boolean networks (AND-nets) on undirected dependency 
graph. That is, on graphs where each edge is bidirectional and hence the 
dependency graph consists of cycles of length two. In particular, their result 
that OR-nets have only fixed points and possibly a limit cycle of length two 
[sil . Lemma 1] follows directly from our results as we explain in Remark 
El 



In [37[, the authors study a smaller family where each edge is undirected 
and each node in the network has a self loop. In particular, they showed 
that disjunctive Boolean networks (which they called OR-PDS) have only 



fixed points as limit cycles [37|, Theorem 3.3]; this follows from our results 
as we show in Remark 12.81 
• Another family consists of the conjunctive Boolean cellular automata and 
was analyzed in In particular, upper bounds for the number of cycles 
are given. Here we study the whole class of conjunctive Boolean networks 
and provide lower and upper bounds for the number of their limit cycles. 
In particular, we present formulas for the exact number of fixed points of 
any conjunctive or disjunctive boolean network, see Equation (|7.4p . 
In this paper we focus only on conjunctive Boolean networks, since for any 
disjunctive Boolean network there exists a conjunctive Boolean network that has 
exactly the same dynamics after relabeling the and 1 and hence we have the 
following theorem. 

Theorem 1.3. Let f = (/i,...,/„) : — > F^' where /, = V • • • V Xi^^ 
be any disjunctive Boolean network. Consider the conjunctive Boolean network 

g = (gi, . . . , Qn) : Fj > F2 , where gi = Xi-^ A - ■ ■ A Xi-, . Then the two phase spaces 

S{f) and S{g) are isomorphic as directed graphs. 

Proof. Let ^ : F'2' — > F^ be defined by ^(xi, . . . , x„) (1 + xi, . . . , 1 + a-„). Then 
it is easy to see that f{xi, . . . ,x„) — [-1 o g o -n)[xi, . . . ,a;„). Thus S{f) and S{g) 
are isomorphic directed graphs. □ 

Remark 1.4. Let G be a graph on n vertices such that the in-degree for each vertex 
is non-zero (G has no source vertex). Then there is one and only one conjunctive 
network / on n nodes such that !D(/) = G. Thus there is a one-to-one correspon- 
dence between the set of all conjunctive Boolean networks on n nodes where none 
of the local functions is constant and the set of directed graphs on n vertices where 
none of the vertices is a source. 

This correspondence was used in [l^ to find the period of a given Boolean 
monomial system (conjunctive Boolean network) and to decide when that system 
is a fixed point system as we will recall in the next section. In this paper we present 
upper and lower bounds on the number of cycles of any length in the phase space 
of any conjunctive Boolean network. Furthermore, we give upper bounds for the 
height. 

In the next section, we recall some results from graph theory as well as results 
about powers of positive matrices that we will use to obtain upper bounds for the 
lengths of transients. 
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2. The Relationship Between Dependency Graph and Dynamics 



Let / : F2 — > F2 be a conjunctive Boolean network, G ~ ^{f) and A the 
adjacency matrix of G. We will assume here and in the remainder of the paper 
that none of the Boolean coordinate functions of / arc constant, that is, all vertices 
of G have positive in-degree. 

2.1. The Adjacency Matrix. Define the following relation on the vertices of G: 
a ~ 6 if and only if there is a directed path from a to 5 and a directed path from 
b to a. It is easy to check that ~ is an equivalence relation. Suppose there arc t 
equivalence classes Vi , . . . , Vj . For each equivalence class Vi , the subgraph Gi with 
the vertex set Vi is called a strongly connected component of G. The graph G is 
called strongly connected if G has a unique strongly connected component. 

There exists a permutation matrix P that permutes the rows and columns of A 
such that 

'Ai Ai2 ■■■ Alt 
A2 ■■■ A2t 



(2.1) 



PAp-' = 











ij represents the 
The form in (|2T]) is 



where Ai is the adjacency matrix of the component Gi, and A, 
edges from the component Gi to Gj, see [SSl, Theorem 3.2.4 
called the Frobenius Normal Form of A. 

Remark 2.1. The effect of the matrix P on the dependency graph is the relabeling 
of the vertices of G such that the diagonal blocks correspond to strongly connected 
components of G. In Example 11.11 above . the adjacency matrix of the dependency 
graph is in the normal form. 



Example 11.11 (Cont.). The dependency graph G in Figure [T] has 3 strongly con- 
nected components and their vertex sets are: Vi = {1,2,3},V2 = {4}, and V3 = 
{5,6}, see Figure [3] (left) . 

For any non-empty strongly connected component Gi, let hi be the conjunctive 
Boolean network with dependency graph 'D{hi) = Gi. Let h : — > Fg be the 
conjunctive Boolean network defined by ft, = (hi, . . . ,ht). That is, the dependency 
graph of h is the disjoint union of the strongly connected graphs Gi, . . . , Gt- 

£^a;amp^e ll.ll (Cont.). The conjunctive Boolean network h corresponding to the dis- 
joint union is ft : F^ — > Fj and given by ft(a;i, . . . , xg) = (fti(a;i, X2, X3), h2{x4), h^{x^, 
where hi{xi,X2,xz) = {x2X3,xi,X2), h2{x/i) = Xi, and h^{x5,X(,) = {xq,X5). 
Now define the following relation among the strongly connected components Gi, . . . ,Gt 
of the dependency graph 23(/) of the network /. 

(2.2) Gi ^ Gj if there is at least one edge from a vertex in Gi to a vertex in Gj. 

Since Gi, . . . ,Gt are the strongly connected components of ©(J), the set of strongly 
components with the relation ^ is a partially ordered set V. In this paper, we relate 
the dynamics of / to the dynamics of its strongly connected components and their 
poset V. 

Example 11.11 (Cont.). The poset of the strongly connected components of / is in 
Figure IHJright). 
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Figure 3. The strongly connected components of / (left) and 
their poset (right). 

For any non- negative matrix A, the sequence {A, A^, . . .} has been studied exten- 
sively, see, for example, 3^ 3§|. Next we use the Boolean operators AND and OR 



to examine the sequence of powers of A and infer results about the conjunctive 
Boolean network that corresponds to an adjacency matrix A. 

2.2. Powers of Boolean Matrices. Let A, B he n x n Boolean matrices (all 
entries are either or 1). Define A(E) B such that (^4 (g) B)ij — Vl-=i(^ifc ^ ^kj), 
where V (resp. A) is the Boolean OR (resp. AND) operator. 

Proposition 2.2. Let A,B be as above and let /, g : — > Fj be the two con- 
junctive Boolean networks that correspond to the adjacency matrices A and B, 
respectively. That is, fi = x1'^X2'^ ■ ■ ■ a;"*" and gi — x'-^^x'^'^ ■ ■ ■ a;^'" /or all i. Then 
the adjacency matrix corresponding to f o g is A® B. 

Proof. It is easy to see that, for all i, 

/i(ffi, • ■ • ,5n) = .9i"---5«" 

^QilfcllH \-a.inbnl ^OiiblnH \-ainbnn 

-^1 ' ' ' -^n 

Since we are working over F2, for all 1 < A; < n, wc have — Xk for all positive 
integers q. Thus 

71 

Xk divides fi{gi, . . . , gn) '^=> ^ ciijbjk > 1 

-^=> a,ijgbj„k = 1 for some jo 
-^=> a.ijg Abj„k = 1 for some jo 

71 

^ V ^ ^j'^ = ^ 

^ (v4®B),fe = l. 
Thus the matrix {A ® B) is the adjacency matrix of / o g. □ 
Throughout this paper, we use A^ to denote A ^ ■ ■ ■ (g) A. 
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Corollary 2.3. Let f be a conjunctive Boolean network and let A be the adjacency 
matrix of its dependency graph £>(/)• Then A* is the adjacency matrix of . 

2.3. The Loop Number. An invariant of a strongly connected graph, called the 
loop number, was defined in [isj . We generalize this definition to any directed 
graph. 

Definition 2.4. The loop number of a strongly connected graph is the greatest 
common divisor of the lengths of its simple (no repeated vertices) directed cycles. 
We define the loop number of a trivial strongly connected graph to be 0. The loop 
number of any directed graph G is the least common multiple of the loop numbers 
of its non-trivial strongly connected components. 

Remark 2.5. Let G be a directed graph and A its adjacency matrix. 

(1) The loop number of G is the same as the index of cyclicity of G as in (4ol | 
and the index of imprimitivity of A as in [ssl : Wi\ . 

(2) The loop number can be computed in polynomial time, see [l^l for an 
algorithm. 

(3) If the loop number of G is 1, the adjacency matrix A of G is called primitive, 
see 

Example 11.11 (Cont.). The loop numbers of Vi,V2,V3 are 1,1,2, respectively. See 
Figure [Jljleft). In particular, the loop number of G is 2. 

Definition 2.6. The exponent of an irreducible matrix A with loop number c is 
the least positive integer k such that — A''. 



The following lemma follows from Proposition 12.21 and the definition above. 

Lemma 2.7. Let f be a conjunctive Boolean network and let A be the adjacency 
matrix of its dependency graph D{f). Suppose the loop number of A is c and its 
exponent is k. Then h(/) = k and p(/) divides c. In particular, C{f)i = for 
every i \ c and hence Equation hl.l]) becomes 

(2.3) C(/) = ^G(/).C„ 

i\c 

Examvle 11.11 (Cont.). The phase space S{f) has 2^ vertices, its cycle structure is 
C(/) = 4Ci + IC2. In particular, the period of / is 2 which is the same as the loop 
number of its dependency graph G. 

Remark 2.8. It is clear that if xi appears in fi for all i, then the loop number of 
'S(f) is one and hence S{f) has only fixed points as limit cycles, this was shown in 
[371 . Theorem 3.3]. Also if each edge in the dependency graph is undirected, then 
the dependency graph is made up of simple cycles of length 2 and hence the loop 
number of S(/) is either two or one. Thus S{f) has only fixed points and possibly 
cycles of length two, which was shown in [s^. Lemma 1]. 



For the case when S(/) is strongly connected p(/) = c as was shown in |18l . 
Theorem 4.13]; in particular, S{f) has a simple cycle of length / if and only if I 
divides c. In general, however, this is not the case. 

Example 2.9. Consider the Boolean network / = {x2,xi,X2X5,x^,X4) ■ 1F| — > F|. 
The graph S)(/) has two strongly connected components with loop numbers 2 and 
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3 respectively, and hence the loop number of T){f) is 6. However, it is easy to see, 
using DVD [29j, that C{f ) = 3Ci + IC2 + 2C3, and hence / has no cycle of length 6. 

The exponent has been studied extensively and upper bounds are known, (40l . 
Theorem 3.11] presents an upper bound for the exponent of any irreducible matrix. 
Using the lemma above, we rewrite this upper bound for the height of any Boolean 
monomial system with a strongly connected dependency graph. 

Theorem 2.10. Let f be a conjunctive Boolean network with strongly connected 
dependency graph D{f), and suppose the loop number ofD{f) is c. Then 

( (n- i)» + i, 1; 

Proof. The proof follows from Proposition 12 . 21 above and [13, Theorem 3.11]. □ 

Furthermore, (40l . Theorem 3.20] also implies an upper bound for the height of 
any conjunctive Boolean network. 

Theorem 2.11. Let f be any conjunctive Boolean network on n nodes. Then 
Hf) < -3n + 2. 



Proof. The proof follows from Proposition 12 . 21 above and [40|, Theorem 3.20]. □ 

We close this section with a classical theorem about positive powers of Boolean 
matrices. This theorem has the remarkable corollary that almost all conjunctive 
Boolean networks have a strongly connected dependency graph and, furthermore, 
have only fixed points as limit cycles. 

Theorem 2.12. f3^ . Theorem 3.5.11] Let N{n) be the number of conjunctive 
Boolean networks on n nodes with strongly connected dependency graph and loop 
number 1. Then 

hm — ~ 1. 

n — >oo 2^ 

2 

In particular, since there are 2" different conjunctive Boolean networks on n nodes, 
as n — > 00, almost all conjunctive Boolean networks on n nodes have a strongly 
connected dependency graph and have only fixed points as limit cycles. 

Although there is a wealth of information about powers of non-negative matrices 
such as the transient length or possible cycle length, very little seems to be known 
about the number of cycles of a given length in the phase space and that is the main 
goal of this paper. In the next section we give a complete answer to this problem 
for conjunctive Boolean networks with strongly connected dependency graph. For 
this class of systems, we find the number of cycles of any possible length in the 
phase space. 

3. Networks with Strongly Connected Dependency Graph 

In this section we give an exact formula for the cycle structure of conjunctive 
Boolean networks with strongly connected dependency graphs. Since "almost all" 
conjunctive Boolean networks have this property by Theorem 12.121 one may con- 
sider this result as giving a complete answer for conjunctive Boolean networks in 
the limit. However, in the next section we will also consider networks with general 
dependency graphs and give upper and lower bounds for the cycle structure. 
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Before deriving the desired state space results we prove some needed facts about 
general finite dynamical systems. 

Lemma 3.1. Let f : X — > X be a finite dynamical system, with X a finite set, 
and let u G X be a periodic point of period t. 

(1) If f^{u) = u, then t divides s. 

(2) The period of f^{u) is t, for any j >1- 

(3) ///*(u) = /-'(u), then t divides s — j. 

Proof. Since t is the period of u and /' (u) = u, by definition s > t. Thus s = qt + r, 
where < r < t. Now u = /"(u) = /''(/'?*(u)) f{u). Since r < t and tis the 
period of u, r = and hence t divides s. That proves ([T]). The proof of ([2]) is 
straightforward, since /*(/^(u)) = /^(/*(u)) = /-'(u). Now wc prove Suppose 
s > j. Since f^u) = /^(u), we have f'-'{fHu)) = r{u) = /^(u). Thus, by g]), 
the period of /^ (u) divides s ~ j. But the period of /^ (u) is i, by ([H). □ 

Lemma 3.2. Let f : X — > X be a finite dynamical system. Then, p(/) = c and 
h(/) = d if and only if c and d are the least positive integers such that /"'^'^(u) ~ 
/™(u) for all m> d and u £ X. 

Proof. Suppose that p(/) — c and h(/) = d. Then for all u G X and m > d, /™(u) 
is a periodic point and hence its period divides c. Thus /'"+^(u) = f'^{f"^(u)) = 
/-(u). 

Now suppose c and d are the least positive integers such that /™+'^(u) = /"(u) 
for all m > d and u £ X. We want to show that p(/) = c and h(/) = d. It is clear 
that /'^(u) is periodic for all u £ X and d is the least such positive number. Thus 
the height of / is d. Also, the period / is c, since c is the least positive integer such 
that /"^(u) = u for any periodic point u £ X. □ 

Theorem 3.3. Let f : X — > X and g : Y — > Y be two finite dynamical systems. 
Define the system h : X x Y — > X x Y by /i(u, u') = (/(u), (jr(u')). Then S{h) = 
Sif) X Sig). 

Proof. This follows from the fact that h{u, u') — (v, v') if and only if /(u) — v and 
.g(u') = v' for all u e X and u' eY. □ 

Corollary 3.4. Let f , g and h be as in Theorem \3.3l Then the period of h is 
p(/i) = lcm{p(/), p((7)} and its height is h(/i) ~ max{h(/), h(g)}. 

Proof. It is easy to see that a set A C X x Y is a cycle in S{h) if and only if 
Ax (resp. Ay) is a cycle in S{f) (resp. S{g)), where Ax .= {u £ X : (u,v) e 
A for some v e Y}. Furthermore, it is clear that the length of the cycle A is the 
least common multiple of the lengths of Ax and Ay. Thus p{h) ~ lcm{p(/), p{g)}. 
The proof of h{h) = max{h(/), h{g)} follows from the definition of height. □ 

Recall Equation in particular, that C{f)m is the number of cycles of length 

TO in the phase space of /, the following corollary follows directly from Theorem 

ESI 
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Corollary 3.5. Let j , g and h he as above. Then, the cycle structure of h is 
C{h) = C{f)C{g) := E,„|p(/i) C{h)mCm, where 

(3.1) C{h)ra= gcd{s,t}C{f),C{g)t. 

t\p{g) 

lcm{s,i}— rn 

That is, the generating function for the cycle structure of h is the product of the 
generating functions of f and g, where Cg ■ Ct = gcd{s,t}C\cm{s.t} ■ 

Let / : — > be a conjunctive Boolean network. Assume that the dependency 
graph of / is strongly connected with loop number c. For any divisor k of c, it 
is well-known that the set of vertices of can be partitioned into c non-empty 
sets Wi, . . . , Wk such that each edge of is an edge from a vertex in Wi to a 

vertex in Wi+i for some i with 1 < i < k and Wk+i = Wi. For a proof of this fact 



see 



38l Lemma 3.4.1(iii)] or [18|, Lemma 4.7]. 
By definition the length of any cycle in the phase space S{f) of / divides c, 
the period of /. For any positive integers p, k that divide c, let A{p) be the set of 
periodic points of period p and let D{k) := A{p). 

p\k 

Lemma 3.6. The cardinality of the set D{k) is \D{k) \ ~ 2^. 
Proof. Let $ ; Ff — > D{k) be defined by 

. ..,Xk) = { xi, . ,xi , X2, . ^. ,X2 , . . 

si times S2 times 

where, without loss of generality. Wi = {vi,i, ■ . ■ ,Vi,si} for all 1 < i < fc. We 
will show that $ is a bijcction. First we show that (f> is well-defined. Let z = 
^{xi, . . . ,Xk), by [13, Theorem 4.10], there exists a positive integer m such that 

ymfc+j(^^-j ^ ( Xj + i, . ^. , Xj+uXj+2, ■ -J , Xj+2j ■ ■ ■ , Xj,. . ^.,Xj ). 

si times S2 times s^ times 

In particular, f^{z) = f^^+h[z) = z. Thus z = ^{xi, . . . ,Xk) e D{k). Since it is 
clear that 4> is one-to one, it is left to show that $ is onto. 

Let z G D{k) which means that f''{z) = z. It is enough to show that Zi^h = Zi^g 
for all 1 < /i, 5 < Si^ and 1 < i < k. Suppose not. Without loss of generality, let 
Zi./i ~ 1 and Zi_g ~ 0. Since k divides c, fiz) = z. But there is a path from vi^h and 
Vi^g of length divides c, contradiction. Hence <!> is a bijection and |£'(fc)| = 2*^. □ 

Corollary 3.7. If p is a prime number and p^ divides c for some k > 1, then 

Proof. It is clear that |-D(1)| = 2. Namely, (0, . . . , 0) and (1, . . . , 1) are the only 
two fixed points. Now if p is prime and A: > 1, then the proof follows from the fact 
that D{p'') = D{p'^^^)[ij A{p''), where 1+J is the disjoint union. □ 



The following theorem gives the exact number of periodic points of any possible 
length. 
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Theorem 3.8. Let f be a conjunctive Boolean network whose dependency graph is 
strongly connected and has loop number c. If c ~ 1, then f has the two fixed points 
(0, 0, . . . , 0) and (1,1,..., 1) and no other limit cycles of any length. If c > 1 and 
m is a divisor of a, then 

1 1 

(3.2) \A{m)\ = ^ ... ^(-l)'i+*2+-+»-2P^""pJ'^"'-f'''"", 

where m = Jl'i=iPi^' prime factorization of m, that is pi, . . . ,pr are distinct 

primes and ki > 1 for all i. 

Proof. The statement for c = 1 is part of the previous corollary. Now suppose that 
c > 1. For 1 < J < r, let m, = pf'' ]T^=l.^^J ■ Then D{m) = A(m) l+J(U;^i 
where 1+J is a disjoint union, in particular, 

r 

A{Tn) D(m) \ |J Dirai). 

The formula l3 . 21 follows from the inclusion-exclusion principle and the disjoint union 
above. □ 

Corollary 3.9. If m divides c, then the number of cycles of length m in the phase 

space of f is C(f)m ~ ^- — ^ — Hence the cycle structure of f is 

m 

c(/) = E^c.- 

^ — ' TO 

m|c 

Remark 3.10. Notice that the cycle structure of / depends on its loop number only. 
In particular, a conjunctive Boolean network with loop number 1 on a strongly con- 
nected dependency graph only has as limit cycles the two fixed points (0,0,..., 0) 
and (1,1,..., 1), regardless of how many vertices its dependency graph has. 

4. Networks with general dependency graph 

Let / : F2 — > F2 be a conjunctive Boolean network with dependency graph 
J)(/). Without loss of generality, the adjacency matrix of T){f) is in the form (j2.1|l . 
Let Gi,...,Gt be the strongly connected components of Ti{f) corresponding to 
the matrices Ai, . . . , At, respectively. Furthermore, suppose that none of the Gi is 
trivial (i.e., Ai is the zero matrix). For 1 < i < t, let hi be the conjunctive Boolean 
network that has Gi as its dependency graph and suppose that the loop number of 
hi is li. In particular, the loop number of / is Z := lcm{/i, . . . , It}. 

In the remainder of the paper, we present lower and upper bounds for the number 
of cycles of a given length. We use the strongly connected components of the 
dependency graph and their poset to infer the cycle structure of /. 

Let Gi and G2 be two strongly connected components in £>(/) and suppose 
Gi ^ G2. Furthermore, assume the vertex set of Gi (resp. G2) is {xi^^, . . . , Xi^} 
(resp. {xj^ ,Xj^}). Without loss of generality, let be a directed edge 

in £>(/) and let S)' be the graph D(/) after deleting the edge (xi^,Xj^). Let g be 
the conjunctive Boolean network such that Tl{g) = D'. 

Theorem 4.1. Any cycle in the phase space of f is a cycle in the phase space of 
g. In particular C{f) <C{g) componentwise. 
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Proof. Let C := {u, /(u), . . . , /™^^(u)} be a cycle of length m in S{f). To show 
that C is a cycle in S(g), it is enough to show that, whenever Xij^-value in u is 0, 
there exists Xj^ g G2, such that there is an edge from Xj^ to Xj-^ and the -value 
in u is 0. Thus, the value of Xj-^ is determined already by the value of Xj^ and the 
edge (xi-ijXjj) does not make a difference here and hence C is a cycle in S{g). 

Suppose the loop number of Gi (resp. G2) is a (resp. b). Now, any path from 
Xi-^ (resp. U i-, ) to itself is of length pa (resp. qb) where q,p > T and T is large 
enough, see [18|, Corollary 4.4]. Thus there is a path from Xi-^ to Xj^ of length qa+\ 
for any q > T. Also, there is a directed path from Xj-^ to Xj^ of length qb — I for 
any q^T . This implies the existence of an edge from Xi^ to Xj^ of length q{a + b) 
for aU q>T. 

Now u = /'"(u) = /'"'^(u), for all fc > 1. Choose q, fc large enough such that 
q{a + b) = km. Then the value of Xj^^ in (u) is equal to zero, since there is 

a path from Xi-^ to Xj^ of length q{a + b) and the value of Xi-^ is zero. Therefore, 
the value of xj^ in u is zero, since u = /™'^(u) = /''(°+''' (u). □ 

Corollary 4.2. Let f and D{f) be as above. For any two strongly connected 
components in S)(/) that are connected, drop all but one of the edges. LetD' be the 
new graph and let g be the conjunctive Boolean network such that D{g) ~ D' . Then 
any cycle in S{f) is a cycle in S{g). In particular C{f) < C{g) componentwise. 

Notice that there are many different possible g one can get and each one of them 
has the same posct as / and provides an upper bound for the cycle structure of /. 
However, we do not have a polynomial formula for the cycle structure of g. When 
we delete all edges between any two strongly connected components we get an easy 
formula for the cycle structure for the corresponding system as we describe below. 
In Section [7] we will present an improved upper bound for the cycle structure of /. 

Let h : F2 — > be the conjunctive Boolean network with the disjoint union 
of Gi, Gt as its dependency graph. That is, h = (hi, . . . ,ht). For h, there are 
no edges between any two strongly connected components of the dependency graph 
of the network. Its cycle structure can be completely determined from the cycles 
structures of the hi alone. 

Theorem 4.3. LetC{hi) = cycle structure of hi. Then the cycle 

structure of h is C{h) = 0^=1 ^(^«) '^'^'^ number of cycles of length m (where 
m\l) in the phase space of h is 

t 

(4.1) c{hu= ^n°^^- 

3i\li i=l 
lcm{ji,...jf}=m 

Proof. This follows directly from CoroUarv 13.51 □ 

Corollary 4.4. Let f and h be as above. The number of cycles of any length in 
the phase space of f is less than or equal to the number of cycles of that length in 
the phase space of h. That is C{f) < C{h) componentwise. 

Example^^{Cont.). Here h = (/ii, /12), where hi{xi,X2) = {x2,xi) and /i2(a;3, ^4, xs) = 
{x5,X3,X4). It is easy to check that C{hi) = 2Ci + IC2, and C(/i2) = 2Ci + 2C3. By 
Theoremltl C{h) = 4Ci + 2C2 + 4C3 + 2Cg. In particular, C(/) < C(/i). 
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5. Bounds on the cycle structure 

Let /, h, Gi, . . . ,Gt be as above and let k be the loop number of d. Further- 
more, let V be the poset of the strongly connected components. Let be the set 
of all maximal antichains m V. For J C [f] = {1, . . . , i}, let .tj := 11 je/ ''^^'^ ^^^^ 
J:=[t]\J. 

Definition 5.1. Let ^ be a subset of the set of limit cycles of h and let s; be the 
number of limit cycles of length i in A. We define and denote ||^|| by 

\\A\\ :=^5.C.. 

i 

Remark 5.2. Using the inclusion-exclusion principle, if A,Ai, ...tAs are subsets of 
the set of limit cycles of h and A = [J^^ Ai, then, 

Mii= E(-i)'"^'iinA-ii- 

jc[s] jeJ 
Definition 5.3. For any subset J C [t], let 

J- := {k : Gj ^ Gk for some j e J}, J- := {k : Gj h Gk for some j G J}, 

:= {k : Gj -< Gk for some j € J}, and := {k : Gj >~ Gk for some j £ J}. 

A limit cycle C in the phase space of / is Jq (resp. Ji) if the Gj component of C is 
(resp. 1) for all j e J. 

Notice that the sets defined above are closely related to upper and lower order 
ideals in posets, see [42|, p. 100]. 

Example \r^{Gont.). Let J {2}. Then = {2,3}, = {1,2}, = 
{3}, and = {1}. 

Definition 5.4. A limit cycle C in S{h) is J— regular if C is both (J^)o and 
for some maximal antichain J in P. A limit cycle is called regular if it is J— regular 
for some maximal antichain J. That is, an element of a J— regular limit cycle 
has in all entries corresponding to strongly connected components lying above 
components in J and 1 in all entries corresponding to components lying below 
components in J. 

Lemma 5.5. Let J be a set of maximal antichains. For J ^ J , let Aj be the 
set of all J— regular limit cycles in the phase space of f . Then, 

II fi Aj\\ = n 

Proof. It is easy to see that a regular limit cycle C G H j^j Aj if and only if the 
only non trivial components of C are in j G H je j ^ 

6. A Lower Bound 

It is easy to see that a regular cycle is actually a limit cycle in the phase space 
of /, and hence the numbers of regular cycles of different length provide lower 
bounds for the number of limit cycles of different lengths in the phase space of /. 
Next we express these numbers in terms of the cycle structure C{hi) of the strongly 
connected components. 
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Lemma 6.1. The cycle structure of the regular limit cycles in S{h) is 

(6.1) E(-i)''"^' n 

Proof. By definition, tlic set of regular limit cycles is U/eo^-'- Now, using the 
inclusion-exclusion principle and Lemma [531 it follows that 

[jAj = Y.{-ir^+'\\f]Aj\\ 

Jen jcn Jej 

□ 

Theorem 6.2. Consider the function 

£(xi,...,xt):= E(-l)''^'^' n 

Then, for any conjunctive Boolean network f with subnetworks hi,. . . ,ht and n 
it's set of maximal antichains in the poset of /, we have 

(6.2) £{C{hi),...,C{ht))<C{f). 

Here, the function C is evaluated using the "multiplication" described in Corollary 
\3.5l This inequality provides a sharp lower bound on the number of limit cycles of 
f of a given length. 

Proof. The inequality follows from the previous lemma. Now to show that this lower 
bound is sharp, it is enough to present a Boolean monomial dynamical system / 
such that C(/) = £(C(/ii), . . . , C{ht)). 

Let / be such that whenever there is a directed edge from Gi to Gj in the poset 
p, there is a directed edge from every vertex in G,; to every vertex in Gj in the 
dependency graph of /. It is enough to check that any cycle in the phase space of 
/ is regular. Let C be a cycle in the phase space of /. Let be the projection of C 
on the strongly connected component Gj. It is clear that is cycle in the C{hj). 
Suppose that is non trivial. Now if Gi < Gj, then all entries corresponding to 
Gi must be one. On the other hand, if Gi ^ Gj, then all entries corresponding to 
Gi must be zero. Therefore, C is regular. □ 

Note that the left side of the inequality (|6.2p is a polynomial function in the 
variables C(hi), with coefficients that are functions of the cycle numbers of various 
lengths of the hi and the anti-chains of the poset of strongly connected components. 
That is, the sharp lower bound is a polynomial function depending exclusively on 
measures of the network topology. 

7. An Upper Bound 

Next we present a polynomial whose coefficients provide upper bounds for the 
number of possible limit cycles in the phase space of /. Similar to the lower bound 
polynomial, this upper bound polynomial is in terms of the strongly connected 
components of / and their poset. However, this upper bound is not sharp. In 
fact, we will show that it is not possible to give a sharp upper bound that is in 
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a polynomial form with constant coefficients. First we prove the following set- 
theoretic equality which we need in this section. 

Lemma 7.1. Let {Ai : i G A} be a finite collection of finite sets of limit cycles in 
S{h), A = flj and B, C A,. Then, for K C A, 

II fl {Au\Bk)\\= J2 (-1)1^1 ||A n nS/.-ll- 

k&K LCK k£L 

Proof. To begin with, observe that 

n {A,\B,)= fl {A\Bu) 

keK keK 

= A\{j (AnSfc). 

kt^K 

Then, using the inclusion-exclusion principle we obtain: 
II fl {Au\Bu)\\ = \\A \ U {Ar^B,)\\ 

keK keK 



LCK keL 

= ^(-i)i^i||Anni3,| 



LCK keL 



□ 



Definition 7.2. A limit cycle C in S{h) is called admissible if, for any j G [t], 

(1) If the coordinates of states in C corresponding to component Gj are 0, then 
C is (j-)o; and 

(2) If the coordinates of states in C corresponding to component Gj are 1, then 
C is 

Notice that any limit cycle in the phase space of / is admissible. In particular, 
we have the following lemma. 

Lemma 7.3. The number of admissible limit cycles in S{h) is an upper bound for 
the number of limit cycles in S{f). 

For J C [t], let Zj be the set of all Jq limit cycles and let Oj be the set of 
all Ji limit cycles. Let K,L C [t], then it is easy to see thatZ^ H Zl — Zkul, 
Ok n Ol = Okul, and the cycle structure of Zk fl Ol is 

\\ZKnOL\\^{K,L) II C{h,), 

where 

'O, if A' nL 7^0, 
1, ifA'nL = 0. 



{K,L) = 



Theorem 7.4. The cycle structure of the admissible limit cycles in S{h) is equal 
^civc[t] keiNUJ^' 

JCMC[t] 

where In = I- U N, j"' ^ U M. 
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Proof. By definition, it follows that the set of admissible limit cycles is 

^dm := {limit cycles in |J ) U (0{j}\0{j}:< ). 

Using the inclusion-exclusion principle, it follows that 
(7.1) 

\\Adm\\ = n E n (%}\%}0n n (0{,}\0{,}.)ii. 

MC[t] 

On the other hand, since HieAr ^{i] — ^n, PljeA/ ^{j} = Om, Zi C .^{^jt and 
Oj C 0{j}r< it follows from Lemma mi that 



= E 

7CA 

= E 

I<ZN 
JCA, 

= E 

ICN 
JCA, 

= E 

7CA 
JCAi 

= E 

7CA 

- E 



-1) 



\{ZNnOM)n{f] Z[,yt n fl o 
\Zn ciOm n Zjy nOj^W 
|Z/„ nOjA/ll 



ICN 
JdM 



Equation ()7.ip then becomes 



\\Adm\\ = II dh,) + y: (-1)1^1+1^^1 e n 



/ce[t] 



JVC[t] 
A7C[i] 



Y[c{hk)+ E (-1) 



JCM 



|Ar| + |A/| + |/| + |,7|/r jM 



{In, J'') n ^('^fc) 



fee[t] 



■fCAfC[t] 
JCMC[t] 
NUM=^0 



= ^ (-i)i^i+i'^^i+i^i+i''i(/^,j^o n c(/».). 



-fC7VC[t] 
JCAfC[t] 



fce/jvUJ*^ 



□ 
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Corollary 7.5. For K,L C [t], let 

0(;,,,)(.xi,...,xO = | ^_ ,;^nL = 0. 

Consider the function 

(7.2) U{x,,...,xt)= 

ICN, JCM 

then 

(7.3) C{f)<U{C{h,),...,C{ht)), 

which provides therefore an upper bound for the cycle structure of f . 

This upper bound is clearly not sharp, since it is easy to find admissible limit 
cycles in the phase space of h that are not cycles in the phase space of /. The 
advantage of this upper bound, however, is that it is a polynomial in terms of the 
strongly connected components and their poset. Furthermore, we will show next 
that this bound gives the exact number of fixed points. 

Theorem 7.6. Let f be a conjunctive Boolean network. Then any admissible fixed 
point is regular and hence the number of fixed points C(/)o in the phase space of f 
is 

(7.4) C(/)o= ^ (-1)1^1+1 2in.e^''l. 

Proof. It follows from the definitions that any admissible fixed point is regular, and 
hence U{2Ci, . . . , 2Ci) = £(2Ci, . . . , 2Ci). In particular, the number of fixed points 
C(/)o is the coefficient of Ci in C{2Ci, . . . , 2Ci). By Equation (|6.ip . we get 

/:(2Ci,...,2Ci) = 5] (-1)1-^1+^ n 2Ci 
= ;^(-i)i-^i+^2in.7e.7-^iCi. 

jcn 

□ 



7.1. An example. Each system g from CoroUarv 14.21 has the same poset as / 
and the cycle structure of each g is an upper bound for the cycle structure of 
/. However, using the poset alone, or using the cycle structure of the strongly 
connected components alone, one can not expect to find a polynomial form with 
constant coefficients that would provide a sharp upper bound, as we now show. 

Let p be the poset on two nodes Gi and G2 where Gi :< G2. Let / be a 
conjunctive Boolean network with a dependency graph that has only two strongly 
connected components which are connected by one edge. Let F be the set of all such 
systems. We will show that there is no polynomial W{xi,X2) = J2i 

n such 

that, for all f & F, W{C{hi),C{h2)) ~ C{f), where hi and /i2 are the conjunctive 

Boolean networks corresponding to the two strongly connected components of /. 

Suppose the loop number of hi is 1 and that of /12 is q where q is prime. Then 

29 — 2 

C(hi) = 2Ci and C(/i9) = 2Ci H C„. It is easy to check that C(f) = 3Ci + 

q 

^—^C„. Since C(f) = W(C(hi),C(h2)), we have 
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2'' — 2 2"^ — 2 

C„ + 3Ci = W{2Ci , 2Ci + C„) 

q q 

= ^a,,(2Cir(2Ci + ^^C,y 

= ^2V„(2Ci + c,y 

913 _ 2 

= [2^Ci + [(2 + y - 2^]C,] 

= J2 2'+^ayCi + ^ 2^a,, [(2 + ?l-^y - 2^]Cg. 

Equating coefficients, for any q, we fiave = X]ij 2'ay [(2 + ^^^)^ - 2-']. 

Therefore, = for all i,j > 1 and hence W{xi,X2) must be of the form 
W{xi,X2) = aiiXiX2 + aiQXi + aoiX2 — aoo- 

Now suppose the loop number of hi is p and that of /12 is q. Then any cycle in 
the phase space of / is regular and hence it is easy to check that C(/) = C{hi) + 
C(/i2) — Ci. In particular, an = and hence W{xi,X2) = aioa^i + aQiX2 — aoo- 

However, for the case when the loop number of hi is 4 and that of /12 is 6, we 
have C(/) = 3Ci + 3C2 + 2C3 + 4C4 + llCg + 2Ci2, C{hi) ^ 2Ci + C2 + 3C4 and 
C(/i2) = 2Ci +C2 + 2C3 + 9C6. In particular, C(/) 7^ VK(C(/ii), C(/i2)). 



8. Discussion 

In this paper we have focused on the class of conjunctive and disjunctive Boolean 
networks and have treated the problem of predicting network dynamics from the 
network topology. Such networks are entirely determined by the topology of their 
dependency graphs, so it should be possible in principle to extract all information 
about the dynamics from this topology. This problem has been the subject of many 
research articles and been solved for the class of XOR Boolean networks, that is, 
networks where all nodes use the XOR Boolean operator. Determining the dynam- 
ics of a network from its topology has practical importance, for instance in the 
study of dynamical processes on social networks, such as the spread of an infectious 
disease. Public health interventions, such as quarantine of selected individuals or 
closure of certain institutions or means of travel, often attempt to alter network 
dynamics by altering network topology. The current study can be seen as a the- 
oretical study that begins to elucidate the role of certain features of the network 
topology in supporting certain types of dynamics. Another example is discussed in 
, where it is shown that the dynamics of a certain gene regulatory network in the 
fruit fly Drosophila melanogaster is determined by the topology of the network. 

We have shown that if the dependency graph of the network is strongly con- 
nected, then the key topological determinant of dynamics is the loop number, or 
index of imprimitivity, of the graph, and one can describe the cycle structure ex- 
actly. The approach here is to look at powers of the incidence matrix of the graph 
in the max-plus algebra. In light of Theorem 12.121 one could argue that this is the 
prevalent case for large networks. What is more, this theorem shows that large 
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networks tend to have loop number one, so that their dynamics is extraordinarily 
simple, consisting only of two fixed points. 

For general dependency graphs, we give lower and upper bounds on the cycle 
structure in the form of a polynomial function determined by structural measures 
of the graph. The lower bound is sharp. 

In some sense, the main contribution of this paper lies in the identification of a 
class of networks that provide an excellent case study of the relationship between 
network topology and dynamics, which can provide hypotheses for more general 
classes of networks. 
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